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We have investigated antidot lattices, which were prepared on exfohated graphene single layers via 
electron-beam lithography and ion etching, by means of scanning Raman spectroscopy The peak 
positions, peak widths and intensities of the characteristic phonon modes of the carbon lattice have 
been studied systematically in a series of samples. In the patterned samples, we found a systematic 
stiffening of the G band mode, accompanied by a line narrowing, while the 2D mode energies are 
found to be linearly correlated with the G mode energies. We interpret this as evidence for p-type 
doping of the nanostructured graphene. 

PACS numbers: 63.20.kd 63.22.Rc 78.30.-j 



Since the first report about the preparation of 
graphene single layers via mechanical exfoliation 
technique^, the interest in this seemingly ideal two- 
dimensional system has grown enormously. For the 
most part, this is motivated by the vision that the two- 
dimensional carbon system is potentially a promising 
candidate for future electronic devices^. Scanning Ra- 
man spectroscopy has proven to be a powerful technique 
for the characterization and investigation of graphene 
samples. It enables a fast and nondestructive investiga- 
tion of the electronic and structural properties of layered 
and laterally-structured samples with sub-micron spatial 
resolution. Via Raman spectroscopy, e.g., one is able to 
identify the number of layers in a sample^'^, to charac- 
terize the edge chirality^ or to probe strain^. Most im- 
portantly for this work, it was shown that the doping in 
graphene can be monitored by Raman spectroscopy^'^, 
even in the case of unintentional doping due to charged 
impurities^—. The most prominent phonon modes in 
two-dimensional carbon lattices are: (i) The G band, 
which stems from inplane LO phonons with E2g sym- 
metry at the center of the Brillouin zone, and which is 
centered around 1580 cm~^ in intrinsic graphene. (ii) 
The D mode (around 1340 cm~^), which is attributed to 
inplane TO phonons from around the K or K' points of 
the Brillouin zone, and which is forbidden due to wave- 
vector conservation. It requires a defect-induced scatter- 
ing process to be observable, and appears in high-quality 
graphene layers at the edges, onlyi^. (iii) The 2D mode, 
which is an overtone of the D mode. Due to the two- 
phonon nature of this mode, it is wave- vector conserv- 
ing. Finally, (iv), the D' band, which is again thought 
to be a defect-induced mode from the maximum of the 
LO phonon dispersion between F and K. In this work, 
we have investigated systematically the Raman modes 
in a series of graphene samples, on which periodic hole 
lattices, i.e., so called antidot lattices, were prepared by 
electron-beam lithography and subsequent ion etching. 
Recently, in transport experiments on similar structures 
two of the authors observed pronounced weak localiza- 
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FIG. 1: (a) Microscope image of a graphene flake. The dashed 
rectangles mark the regions, where antidot lattices were pre- 
pared, (b) Scanning electron micrograph of region B, and, (c) 
of region A. (c) Representative Raman spectra of regions A, 
B, and R on the sample. 



tion due to strong intervalley scattering at the antidot 
edges and found a transport gap by analyzing thermally 
activated transport around the Dirac point. Here, we 
find in our Raman experiments a strong correlation of 
the energetic positions of the G modes and 2D modes 
in all antidot samples, from which we conclude that the 
nanoscale patterning produces a p-type doping. We have 
found doping levels in the range of 3 — 7 x 10^^ cm~^. 

The graphene samples were prepared on Si wafers with 
300 nm Si02 cap layers by the well-known mechanical 
exfoliation technique^. On selected single layer flakes, 
antidot lattices with periods of 80 nm, 100 nm, 200 nm, 
and 400 nm were prepared using electron beam lithogra- 
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phy and oxygen reactive ion etching. The hole diameters 
are between about 50 nm and 60 nm. In this series of 
samples, we have intentionally deposited no metal con- 
tacts in order to avoid unwanted additional doping due to 
the presence of metals on the samples^. We use below, 
however, measurements on a gated graphene reference 
sample for estimates of the hole densities. The Raman 
experiments were performed at room temperature with 
a 532 nm diode-pumped solid state laser and a micro- 
scope setup. The samples were mounted on a motor- 
ized x-y translation stage with minimal step sizes of 100 
nm, the laser-spot diameter was 800 nm. Raman spec- 
tra were recorded with a TriVista triple Raman system, 
equipped with a nitrogen-cooled CCD detector. Figure 
[TJi shows an optical microscope image of one of the in- 
vestigated samples. In the regions A and B, indicated by 
dashed rectangles in Fig. [TJi, antidot lattices with a pe- 
riod of 80 nm were prepared. The holes in region A have 
a slightly larger diameter (~ 60 nm, see Fig. [it) than 
the holes in region B (~ 50 nm, see Fig. [Ud). Region 
R is an un-patterned part of the sample. Note that on 
the left edge of region A and on the upper right edge of 
the reference region R, there are multi-layer areas, visi- 
ble by the much darker contrast in the microscope image 
(Fig. [TJi). Typical Raman spectra of the three regions, 
A, B and R, are displayed in Fig. [l]i. In the antidot 
regions, A and B, pronounced D peaks are observable, 
along with D' peaks. We interpret this to be dominant ly 
caused by the introduction of additional edges into the 
sample due to the creation of the holes, and not by ad- 
ditional defects or disorder, which might be caused by 
the preparation process. The reason for this interpreta- 
tion is evident from the scanning Raman images, shown 
in Fig. [2j Figure 2 shows false color plots of data, ob- 
tained from scanning Raman experiments on the same 
sample, as shown in Fig. [H In these measurements, in- 
dividual Raman spectra were taken with lateral spatial 
resolution of about 800 nm, in steps of 500 nm in x and y 
directions. From the raw spectra, peak positions, intensi- 
ties and linewidths of the phonon modes were extracted. 
The image of the G mode intensity (Fig. [2^) essentially 
reproduces the topology of the whole flake. The inten- 
sity of the G mode increases from region A, which con- 
tains the larger holes, to region B, with smaller holes, 
to the un-patterned reference region R. A maximum in- 
tensity is visible in the multi-layer part, left of region 
A. On the other hand, a pronounced intensity of the D 
mode is observed in the antidot areas A and B, only, with 
higher intensities in A, which contains the larger holes, 
i.e., with longer edges (Fig. [2]3). Interestingly, the ener- 
getic position of the G mode is systematically higher in 
the antidot lattices than in the reference part (see Fig. 

While the G mode energy in the reference part is 
between about 1589 and 1590 cm~^, it is in regions A 
and B upshifted to about 1592 to 1594 cm~^. Evidence 
for the high crystalline quality of the prepared samples 
can be deduced from the image of the linewidths of the G 
mode, displayed in Fig. [2li. The linewidths are over the 
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FIG. 2: False color plots of scanning Raman experiments on 
the same sample as in Fig. 1. Displayed are measured images 
of (a) the intensity of the G mode, (b) the intensity of the D 
mode, (c) the energetic position of the G mode, and, (d) the 
line width of the G mode. The numbers at the scale bars in 
(c) and (d) are given in cm""*^. 



whole single-layer region between about 9 and 11 cm~^, 
comparable to pristine graphene samples^^^^, and to plain 
gated graphene outside the Landau-damping region (see 
below Fig. [3b). For significantly disordered samples one 
would expect hnewidths in the range of several tens of 
wavenumbers^^. This clearly excludes the interpretation 
that the observed stiffening of the G mode in the antidot 
sample is caused by disorder. The narrow and almost 
unchanged hnewidths indicate that doping is the reason 
for the upshift, caused by the nonadiabatic removal of 
the Kohn anomaly (see below) The most important 
results of the present work are shown in Fig. [3l Figure 
[3K displays a summary of the most relevant data from 
the whole series of investigated antidot samples. If we 
plot, sample by sample, the positions of the observed 2D 
modes versus the positions of the G modes, we find an ap- 
proximately linear correlation (Fig. [3^^). The linewidths 
of the observed G modes in the different samples tenta- 
tively decrease with increasing G mode energy (inset of 
Fig.[3K). The correlation between G mode and 2D mode 
positions was recently investigated systematically on a 
large number of graphene samples, both, in as-deposited 
and in gated sample^i^^. It was reported that while the 
2D mode position decreases with increasing G mode po- 
sition in n-doped samples, it increases with increasing G 
mode position in p-doped graphene^^. The dashed line 
in Fig. [3^1 indicates, for the range of G mode energies rel- 
evant for our work, the slope of the approximately linear 
correlation, as reported in Ref^^ for p-doped graphene. 
In stark contrast, for n-type doping the slope should be 
negative. One can see a fairly good agreement of the 
data, extracted from our measurements on the antidot 
samples, with the dashed line in Fig. [3^. Please note 
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FIG. 3: (a) Plot of the energetic positions of the observed 
2D modes versus the positions of the G modes of aU inves- 
tigated antidot samples. The legend indicates the periods of 
the investigated antidot lattices. The diameter of the holes 
is approximately 60 nm. The label (s) means slightly smaller 
holes (~ 50 nm). The inset shows the linewidths of the G 
modes versus energetic position, (b) Plot of the energetic po- 
sitions of the G mode and its linewidth versus gate voltage, 
as measured in a plain, gated graphene sample. 



that the 2D mode energies in our experiments are slightly 
lower than in Ref.^^, since we measured with a 532 nm 
line (514 nm line in Ref.^^): The 2D mode energy in- 
creases with decreasing laser wavelength due to the dou- 
bly resonant intervalley scattering process^^^. Therefore, 
the dashed line in Fig. [3^1 is rigidly shifted by about 17 
cm~^ to lower energies as compared to 2D mode ener- 
gies m Fig. 3 of Ref.^^ Obviously, our data in Fig. [3^ 
suggest a p-type doping of the antidot lattices. We note 
here, that in previous Raman investigations of graphene 
edges, also indications of a stronger doping in close vicin- 
ity to the edges, as compared to the bulk of the samples, 
were found^^. Based on this, and on our observations, 
we suspect that the doping effect reaches up to a max- 
imum distance of a few tens of nanometers around the 
edges of the holes. Hence, the amount of p doping ap- 
pears to increase systematically with the ratio of etched 
to unetched graphene area, which is tentatively verified 



by the data in Fig. [3^1. As a result, the plain reference 
sample regions are not affected by the doping. In order 
to quantify the amount of p doping, we have compared 
our results to measurements on an un-patterned gated 
sample, where doping can be introduced electrostatically 
(cf. Refs.^^^). Figure [3]3 shows the observed G mode posi- 
tions and linewidths versus applied gate voltage for elec- 
trostatic p doping. The gate voltage where the Fermi 
energy is at the Dirac point is indicated by a vertical 
blue dashed line. The G mode is Landau-damped due to 
decay into vertical electron- hole excitations, and hence 
has a larger linewidth, as long as the Fermi energy is 
smaller than half of the G mode energy (gray-shaded re- 
gion in Fig.[3]3)^. In accordance with Ref^^, the observed 
G mode energy in the gated sample increases almost lin- 
early with increasing negative gate voltage (indicated by 
dashed black line in Fig. [SJd) due to the nonadiabatic 
removal of the Kohn anomaly at the T point In the 
low-density range, within the Landau-damped region, the 
disorder potential dominates the charge state^^, leaving a 
nonzero average doping, and, hence, an almost constant 
G mode energy. From transport experiments on similar 
samples we know that the relation between hole (elec- 
tron) density and applied gate voltage, measured with 
respect to the Dirac point, is p ^ 7.2 x 10^^ cm~^x Vgate 
[V] . From this and the measured slope of G mode energy 
versus gate voltage in Fig. [3]3, we conclude that the G 
mode energy increases by 1 cm~^ in the graphene samples 
per ~ 4.5x10^^ cm~^ increase in hole density, starting at 
a G mode energy of 1580 cm~^ at zero density. Of course, 
this relation holds for densities > 1 x 10^^ cm~^, above 
the disorder-dominated regime, only. With this calibra- 
tion we can assign hole densities between about 3 x lO"*^^ 
cm~^ (400 nm sample in Fig. [3ti) and 7 x 10^^ cm~^ (80 
nm sample in Fig.[3ti) to our investigated antidot samples 
in Fig. E^. 

In conclusion, we have investigated antidot graphene 
samples by scanning Raman spectroscopy. The overall 
crystalline quality of the antidot lattices is as good as 
that of the as-deposited layers. An observed systematic 
stiffening of the G mode in the antidot samples, accompa- 
nied by a similar increase in 2D mode energy lead to the 
conclusion that the introduction of the antidot lattices 
causes a p-type doping. In all likelihood, the observed p 
doping is a consequence of the patterning process here, 
and not a result of the specific antidot pattern. 
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